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ABSTRACT 

The buckling strength of thin-walled 
stiffened or unstiffened, metallic or composite 
shells is of major concern in aeronautical and 
space applications. The difficulty to predict the 
behavior of axially compressed thin-walled 
cylindrical shells continues to worry design 
engineers as we enter the third millennium. 
Thanks to extensive research programs in the 
late sixties and early seventies and the 
contributions of many eminent scientists, it is 
known that buckling strength calculations are 
affected by the uncertainties in the definition of 
the parameters of the problem such as 
definition of loads, material properties, 
geometric variables, edge support conditions 
and the accuracy of the engineering models 
and analysis tools used in the design phase. 

The NASA design criteria monographs [1] 
from the late sixties account for these design 
uncertainties by the use of a lump sum safety 
factor. This so-called “empirical knockdown 
factor y” usually results in overly conservative 
design. Recently new reliability based 
probabilistic design procedure for buckling 
critical imperfect shells have been proposed 
[2,3]. It essentially consists of a stochastic 
approach which introduces an improved 
“scientific knockdown factor X a '\ that is not as 
conservative as the traditional empirical one. 

In order to incorporate probabilistic 
methods into a High Fidelity Analysis Approach 
one must be able to assess the accuracy of the 
various steps that must be executed to 
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complete a reliability calculation. In the present 
paper the effect of the size of the experimental 
input sample on the predicted value of the 
scientific knockdown factor X a calculated by 
the First-Order, Second-Moment Method is 
investigated. 

INTRODUCTION 

Buckling strength of thin-walled stiffened 
or unstiffened, metallic or composite shells is 
of major concern in many aeronautical and 
space applications. It is well known that the 
critical buckling load is affected by the 
uncertainties in the definition of loads, material 
properties, geometric variables, engineering 
models and the accuracy of the analysis tools 
used in the design phase. The NASA design 
criteria from the late sixties account for these 
design uncertainties by the use of a lump sum 
safety factor, the so-called “knockdown” factor 
y[1], which usually results in an overly 
conservative design. 

Lately, probabilistic design procedures 
have been proposed as a viable alternative 
[2,3]. It is felt that quantifying and 
understanding the “problem uncertainties” and 
their influence on the design variables provides 
an approach which will ultimately lead to a 
better engineered, better designed and safer 
structures. 

However, the vast majority of practicing 
engineers agree that true reliability must be 
demonstrated and not just estimated from 
analysis. It is the authors’ opinion, that before 
the engineering community will begin in large 
numbers to accept the current generation of 
probabilistic tools, two conditions must be 
satisfied. First, there must be test-constructed 
data bases which can help in mapping the 
input parameter uncertainties into probability 
density functions. In addition, there must be 
failure data bases, which can be used for test 
verification of the probabilistic failure load 
predictions. 
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It is generally agreed that, in order to make 
the development of the Advanced Space 
Transportation System a success and to 
achieve the very ambitious performance goals 
(like every generation of vehicles lOx safer and 
lOx cheaper than the previous one), one must 
make full and efficient use of the technical 
expertise accumulated in the past 50 years or 
so, and combine it with the tremendous 
computational power now available. It is 
obvious that with the strict weight constraints 
used in space applications these performance 
goals can only be achieved with an approach 
often called “high fidelity analysis”, where the 
uncertainties involved in a design are 
simulated by refined and accurate numerical 
models. In the end the use of “high fidelity” 
numerical simulation will also lead to overall 
cost reduction, since the analysis and design 
phase will be completed faster and only the 
reliability of the final configuration needs to be 
verified by structural testing. 

In order to incorporate probabilistic 
methods into a High-Fidelity Analysis 
Approach one must be able to assess the 
accuracy of the various steps that must be 
executed to complete a reliability calculation. 

The central problem in the application of 
stochastic processes is the estimation of the 
various statistical parameters in terms of real 
data. According to Papoulis [4], using 
ensemble averaging to evaluate the lower 
order statistical moments is appropriate if a 
sufficiently large number of realizations of the 
random vector X(m) 

are available. 

This paper deals specifically with this 
problem. To investigate the effect of the size of 
the input sample on the predicted value of the 
scientific knockdown factor X a , the First- 
Order, Second-Moment Method [3] is applied 
successively to sample groups of different 
sizes. 


TEST PROGRAM 

Using STONIVOKS [5], a special purpose 
testing machine developed at the Structures 
Laboratory of the Faculty of Aerospace 
Engineering in Delft, a sample of 32 nominally 
identical seamless stainless steel “beer cans” 
were tested. Figure 1 shows a typical test 
specimen before top and bottom are cut off. 
Figure 2 gives the general dimensions of a 
typical specimen and shows the variation of 
the wall thickness along 4 equally spaced 
generators. Notice the increasing wall 
thickness towards the open end of the can. 


The imperfection surveys of the test- 
specimen were carried out with the 

STONIVOKS testing machine of the Delft 
University of Technology (see Fig. 3). The 
testing machine mainly consists of a rotating 
platform on which the test specimen is placed, 
and a vertical moving pick-up (LVDT 
transducer). The shell is mounted between two 
circular end discs. The top and bottom ends of 
the can are cut off giving a cylinder of length 
100 mm. Next this specimen is placed between 
the end disks in a circular channel which is 
filled with molten “Cerrobend”. When the 
“Cerrobend” solidifies the edges of the cylinder 
are fully clamped (Fig. 4). The rotary moment 
of the platform and the vertical movement of 
the pick-up are synchronized in such way that 
one revolution of the specimen corresponds to 
a vertical displacement of the pick-up by 1 mm. 
The number of measurements in the 
circumferential direction is fixed at 100. Since 
the usefull length of the cylinders is about 80 
mm the number of measurements per test is 
8000. As the rotary movement of the platform 
(and subsequently the vertical movement of 
the displacement pick-up) is continuous, the 
measuring pattern is a helix over the outside 
surface of the specimen. Figure 3 shows the 
specimen installed in its testing position. For a 
detailed discussion of this apparatus see Ref. 
[5]. All buckling tests were preceded by a 
complete imperfection survey controlled and 
recorded by means of a Hewlett-Packard 
HP9825S desktop calculator. 

Data Reduction 

The data reduction process involves 4 
steps, namely: interpolation of the experiment- 
al data, elimination of the rigid body motions, a 
best fit correction, and finally a harmonic 
analysis. 

As stated before, the measuring pattern is a 
helix over the outside surface of the test 
specimen. In order to make the measurements 
suitable for a harmonic analysis the imper- 
fection data must be interpolated in axial 
direction. Considering the error level present in 
the measured data (due to the measuring 
system) a linear interpolation is sufficient. 

Although production and assembly of the 
testing machine and the preparation of the test 
specimen was very accurate, a certain amount 
of rigid body motion of the test specimen was 
unavoidable. This rigid body motion is mainly 
caused by the fact that the center line of the 
specimen never exactly coincides with the axis 
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of rotation of the platform. Because the 
absolute values of the imperfections are small, 
it is essential to correct the measured 
imperfection pattern for these rigid body 
motions. To measure this rigid body motion, a 
pair of transducers (LVDT’s) are used which 
measure the displacements of the outer rim of 
the tyvo end disks. The outer rim of each disk is 
considered to be concentric to a high degree. 
An analysis of this rigid motion shows that the 
displacements measured at the outer rim of the 
end disks may be considered as to be 
sinusoidal. After calculation of the average 
displacements of each rim and subtracting this 
value from the measured displacements on the 
outer rim of the upper and lower ring a linear 
interpolation of the measured imperfections in 
the axial direction is carried out. 

Before calculating the coefficients of the 
double Fourier series representations of the 
measured contours, it is necessary to 
determine what is to be considered the perfect 
shell. This is done by fitting a best-fit cylinder 
(Fig. 5) to the measured data of the initial 
imperfection scan. The method of least 
squares is used to calculate the eccentricities 
Xi, Yi, the rigid body rotations ei, £2 an ^ the 
mean radius R. Finally, the initial imperfections 
are defined by recomputing the measured 
distances with respect to the newly found 
“perfect” cylinder. The recalculated radial initial 
imperfections of the isotropic shell I W 1-1 6 are 
shown in Fig. 6. The coefficients of the 
following double Fourier series, referenced to 
the “best-fit cylinder” 
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are calculated numerically. 

By assuming that the Fourier coefficients of 
the initial imperfections are the basic random 
variables Xj , we have a sample of 32 random 
fields representing the 32 nominally identical 
shells. For further details of the test-data the 
interested reader should consult Ref. [6], 


g(X) = A S (X)-A. 

where A, = P/P C is the suitably normalized 
load parameter, A s is the random collapse 
load of the shell, and the vector X represents 
the Fourier coefficients of the initial 
imperfections. Clearly the response function 
g(X) = 0 separates the variable space into a 
“safe region" where g(X) > 0, and a “failure 
region” where g(X) < 0 . The reliability R(k ) , or 
the probability of failure Pf (^) can then be 


calculated as 

R(^) = 1-P f (A.) (3) 

where 

Pf (A,) = Prob {g(X) < 0} = [ ff x (x)dx (4) 


' g(X)<o 

The limit state function g(X), if so desired, can 
be determined with great accuracy with 
currently available nonlinear finite element 
codes such as STAGS [7]. Flowever, the 
evaluation of the multi-dimensional probability 
integral, where the domain of integration 
depends on the shape of the response (or limit 
state) function, is by no means trivial. 

Using as an approximation the First-Order, 
Second-Moment Method to evaluate this 
integral involves linearization of the response 
function g(X) -» Z(X) at the mean point and 
knowledge of the distribution of the random 
vector X. To combine the use of numerical 
codes with the mean value First-Order, 
Second-Moment Method, one needs to know 
the lower order probability characteristics of Z. 
In the first approximation the mean value of Z 
is 


E(Z) = E(A S ) - X= EmXf,.... X n )] - X (5) 
~'P[E(X 1 ) E(X n )]-A 


whereas the variance of Z is approximated by 
var(Z) = (6) 



STOCHASTIC STABILITY ANALYSIS 

The collapse problem of axially 
compressed isotropic cylinders can best be 
formulated in terms of a response (or limit 
state) function 


where cov(Xj,Xj) is the variance-covariance 
matrix. The calculation of the value of 
\j/[E(Xi),...,E(X n )] , which is the deterministic 
collapse load of the imperfect shell with mean 
imperfection amplitudes, and of the derivatives 


dxy/d^j (or 3v) i/d^) is done with MIUTAM [8], 
the code that was chosen for the numerical 
work. The calculation of the derivatives 
5\|//3^j (or 9v|f/a^|<) is performed 
numerically by using the following numerical 
differentiation formula evaluated at values of 
5j=E(Xj) (or =E(X k )) 

iV = (7 ) 

Ai >\ 


The mean values and the variance- 
covariance matrix of the basic random 
variables Xj are obtained by evaluating the 
following ensemble averages for a sample of 
experimentally measured initial imperfections 

M , M 
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where M is the number of sample shells, and 
m is the serial number of shells. See 
Reference [9] for further details. 

Having obtained the quantities E(Z) = 
E(A S )-A. and Var(Z), one can proceed to 
estimate the probability of failure Pf(A) as 

0 

Pf (A.) = Prob (Z < 0) = F Z (0) = j f z (t)dt (11) 


where a = E(Z) and oz = ^/Var(Z) . Further 


u 

Fz(0)= | fz (t) dt = -^ + erf 
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where p = a/crz is the reliability index and 
4>(P) is the standard normal probability 
distribution function. The error function erf(p) 
is defined as 



Finally, the reliability R(A,) will be estimated as 
R(A,) = 1 - P f (A.) = 1 - Prob (Z < 0) (15) 


= 1 - Fz (0) = ^ + e rf (P) = 4>(P) 


NUMERICAL RESULTS 


For the statistical calculations the data 
associated with the seamless stainless steel 
“beer scans", tested in Delft in 1987 during the 
test program described earlier, are used. The 
shell properties and test results are given in 
Table 1. 

Relying on the results of earlier 
investigations of the buckling behavior of the 
imperfect isotropic shell A-8 [3], it was decided 
to use the following modified 3-mode 
imperfection model of Koiter [10] for the 
collapse load calculations. 

W 7 .x 

— = 5l C0Sl c fJt- + 

L (16) 

(^2 sink-|rt^- + ^3 sink 2 n^)cos^ 


where Fz(t) is the probability distribution 
function and fz(t) is the probability density 
function of Z. 


Assuming that the limit state function Z is 
normally distributed, then 


fz(t) = 



1 t-a 

exp [ ( 

2 o Z 
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where 


^1 - A a xi I ^2 - A asy 



and k^j and k2 are roots of the equation which 
defines the so-called Koiter circle [1 1] 
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The Koiter circle is the locus of a family of 
modes belonging to the lowest eigenvalue 
X c = 1.0 , where by definition 


N, 


N c ^ 


;N C ^ = 


Et 

cR 
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and No is the applied axial compressive stress 
resultant. The circumferential wave number f. 
is chosen such that k-j = 1 . This yields the first 
Koiter triad [8]. The ratios of the amplitudes of 
the imperfections are chosen such that the 
direction of the postbuckling path coincides 
with the path of steepest descent, which yields 
the most adverse imperfection shape [10]. 

If the measured initial imperfection is 
represented by the double Fourier series of Eq 
(1) then its root-mean-square value is by 
definition 


A rms 
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The derivatives dy / are calculated as 

follows. For the increment of the random 
variable in Eq. (7), 1% of the original mean 
value of the corresponding equivalent Fourier 
coefficient is used, so that A^j = 0.01 ■ E(Xj) . 

The calculated derivatives are listed in Table 4. 
In this study the increments of the path 
parameter are chosen in such a way that the 
limit loads are found accurate to within 0.01%. 
Next, using the sample variance-covariance 
matrix displayed in Table 5, one can evaluate 
the mathematical expectation and the variance 
of Z. The results of these calculations are 
E(Z) = 0.550629 - X and Var(Z) = 0. 1 51 08 1 0.~ 3 . 
Finally, the reliability is calculated directly from 
Eq. (15) and is plotted in Fig. 7. Notice that for 
a reliability of 0.99999, one obtains a scientific 
“knockdown” factor X a = 0.49 . Next the same 
calculations must be repeated for the other 
groups. The results of the calculations for the 8 
groups of 4 shells each are summarized in 
Table 6. The results of similar calculations for 4 
groups of 8 shells each and for 2 groups of 16 
shells each are shown in Tables 7 and 8 
respectively. Using all 32 shells as one sample 
group yields the following results: 

E(Z) - 0.552342 - X ; a z = 0.0234049 

and X a = 0.44 for a specified reliability of 
R(?0 = 0.99999 . 


Thus 
(Arms ) 2 



n i n 2 
k=1, f=2 




(20) 


Using the imperfection model of Eq. (16) the 
corresponding equivalent imperfection 
amplitudes are listed in Table 2. 

In order to investigate the effect of the size 
of the input sample on the predicted value of 
the scientific knockdown factor X a initially the 
32 shells are divided in 8 groups of 4 shells 
each. Next for every group a separate 
reliability calculation is carried out using the 
First-Order, Second-Moment Method (FOSM). 
When applying this method to the first group of 
4 shells each the mean buckling load has to be 
calculated first. Using the imperfection model 
of Eq. (16) with the mean values of the 
corresponding equivalent imperfection 
amplitudes listed in Table 3, the result of the 
calculation is E(A S ) = 0.550629 . 


CONCLUSIONS 


A comparison of the “scientific knockdown 
factors X a " obtained with the different sample 
sizes listed in Table 8 indicates, that evaluating 
the mean values and the variance-covariance 
matrices via ensemble averaging of 
experimental data yields accurate results if the 
sample size is 16 or greater. 

Comparing the probabilistic buckling load 


P a = X a P cf = 0.44(-7909.734) = -3480.283 lbs 


which has a probability of failure of 


P f =1-R(>t) = 1-0.99999 = 0.1 10“ 5 

with the experimental buckling loads shown in 
Table 1, it is seen that 3 of the test shells 
buckled below this value. Thus it is evident that 
the simple initial imperfection model of Eq. (3) 
does not model the collapse behavior of the 
shells tested accurately enough. In view of this 


one must conclude that in order to improve the 
accuracy of the probabilistic predictions to a 
level required by a “High Fidelity Analysis”, one 
must investigate the role of using more refined 
mechanical models to calculate the collapse 
loads and employ more advance methods for 
evaluating the multidimensional probability 
integral of Eq. (4). 
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Table 1 Experimental buckling loads and geometric and material properties of the IW1 -shells [6] 



Pexp 


Pexp 


Pexp 


Pexp 

IW1-16 

3.05 

IW1-24 

4.27 

IW1-33 

4.03 

IW 1 -42 

3.82 

-17 

3.53 

-26 

3.99 

-34 

4.68 

-43 

3.83 

-18 

4.50 

-27 

4.16 

-36 

4.43 

-44 

4.23 

-19 

4.51 

-28 

4.24 

-37 

3.55 

-45 

3.99 

-20 

3.89 

-29 

4.49 

-38 

4.20 

-46 

3.35 

-21 

4.01 

-30 

4.46 

-39 

4.00 

-47 

3.51 

-22 

3.82 

-31 

4.47 

-40 

4.08 

-48 

3.43 

-23 

4.50 

-32 

4.01 

-41 

4.03 

-49 

3.48 

For all shells: 

R = 33.0 mm ; t = 

0.1mm ; 

E = 2.08 10 5 

N/mm 2 I 

Pexp in kN 



L = 100.0mm 


v = 0.3 



Table 6 Reliability calculation for groups of 4 shells each 


Group 

E(Z) 


X a 

1 

0.550629-^ 

0.0122915 

0.49 

2 

0.543369-iV 

0.0156275 

0.47 

3 

0.553647-?i 

0.0289698 

0.42 

4 ' 

0.5601 22-?i 

0.0369662 

0.39 

5 

0.540552-?. 

0.0282565 

0.41 

6 

0.536736-A, 

0.0235181 

0.43 

7 

0.571050-?. 

0.0080433 

0.53 

8 

0.566009-?, 

0.0108056 

0.51 


Table 7 Reliability calculations for groups of 8 shells each 


Group 

E(Z) 

°Z 

x a 

1 

0.546944-?. 

0.0136657 

0.48 

2 

0.556835-?. 

0.0307617 

0.42 

3 

0.538439-?. 

0.0241657 

0.43 

4 

0.568504-?. 

0.0092468 

0.52 


Table 8 Reliability calculations for groups of 1 6 shells each 


Group 

E(Z) 

°Z 


1 

0.551 832-?i 

0.0232321 

0.44 

2 

0.552866-?. 

0.0243072 

0.44 


Table 9 Influence of Sample Size on Reliability Predictions 


Groups of 4 shells 

Groups of 8 shells 

Group of 16 shells 

Group of 32 shells 

#1 

#2 

X a =0.49 
X a =0.47 

#1 

?. a =0.48 



#3 

#4 

?. a =0.42 
?. a =0.39 

#2 

?. a =0.42 

#1 7 a =0.44 


#5 

#6 

?i a =0.41 
X a =0.43 

#3 

X a =0.43 



#7 

#8 

?. a =0.53 
X a =0.51 

#4 

X a =0.52 

#2 X a = 0.44 

#1 ?i a =0.44 







Fig. 3 Displacement pick-ups for contour and rigid body motion measurements [5] 



Fig. 4 Test specimen with end disks 


in 



First Order Second Moment Reliability Plot for Group of 4 beer cans 
Wbor/h - Koiter's 3 — modes imperfection model MODIFIED 



La m b d a 

Fig. 7 Reliability of the first group of 4 IWI-shells using 
Koiter’s modified 3-mode imperfection model 


